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Screening is one of the most important concepts in the study of charged systems. Near a dielectric 
interface, the ion distribution in a salt solution can be highly nonuniform. Here, we develop a theory 
that self-consistently treats the inhomogeneous screening effects. At higher concentrations when the 
bulk Debye screening length is comparable to the Bjerrum length, the double layer structure and 
interfacial properties are significantly affected by the inhomogeneous screening. In particular, the 
depletion zone is considerably wider than that predicted by the bulk screening approximation or the 
WKB approximation. For asymmetric salts, the inhomogeneous screening leads to enhanced charge 
separation and surface potential. 
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Screening due to the ionic atmosphere, introduced by 
Debye and Hiickel more than 90 years ago[i|, is one of 
the most important concepts in the study of charged sys¬ 
tems. Screening has profound effects on essentially all 
properties of biophysical and salt-containing soft matter 
systems [3-0]- In a homogeneous bulk solution, screening 
is most commonly manifested as an exponential damp¬ 
ing of the long-range Coulomb interactions between two 
test charges. When the ion distribution is nonuniform, 
as in the vicinity of a charged surface or an interface with 
dielectric discontinuity, screening also becomes inhomo¬ 
geneous. However, in spite of the ubiquity of systems 
with nonuniform ion distributions, a rigorous treatment 
of screening in such systems is still lacking. 

For a salt solution near a dielectric interface, e.g., the 
water/air interface, the repulsive image force creates a 
depletion layer, whose theoretical treatment was pio¬ 
neered by Onsager and Samaras (OS) @ . This problem is 
related to a number of phenomena, such as conductivity 
in artificial and biological ion-channels [TH-QK stability 
of colloidal, bubble and protein suspensions |l5l - [l8j | . and 
the rate of ozone consumption [l^. [2ril |. Assuming that 
the image force is screened by the bulk screening length, 
the OS theory qualitatively explains the excess surface 
tension and yields agreement with experimental data at 
low salt concentrations (c;, < O.OIM). However, there is 
large discrepancy between the OS theory and experiment 
data at high salt concentrations (c^ > 0.1M)[^,[^|. The 
OS theory predicts an ever decreasing width of the de¬ 
pletion layer with salt concentration, which results in a 
concave downwards curve for the surface tension vs. Cb- 
In contrast, experimental data show essentially a linear 
increase of the surface tension. To reconcile this dis¬ 
crepancy, an exclusion zone of constant width or large 
hydration radius of the ion is usually invoked [2^. 

An obvious effect missing in the OS theory[^ and in 
subsequent modifications [231-13J is the spatially varying 
screening of the image force near a dielectric interface (see 
Fig. 1): the ion concentration changes gradually from 
zero at the interface to the bulk value. In this depletion 



FIG. 1: Schematic of the inhomogeneous screening near the 
dielectric interface. The red and yellow spheres represent the 
cations and anions, respectively. The two test ions are la¬ 
beled, with Ion 1 located very close to the interface and Ion 
2 approaching the bulk solution. The two black spheres are 
the image charges corresponding to the two test ions. 


layer, the ionic cloud is highly anisotropic, giving rise to 
different features of the screening near the interface from 
the homogeneous and isotropic bulk. Close to the inter¬ 
face, ions are strongly depleted; the local ionic strength 
around the test ion (Ion 1 in Fig. I) is much lower than 
the bulk. The bulk screening approximation clearly over¬ 
estimates the screening effect and hence underestimates 
the image force. Even for an ion approaching the bulk 
(Ion 2 in Fig. 1), the screening is still weaker than in the 
bulk due to the long-range, accumulative effects from the 
depletion zone. This feature extends the effective range 
of the image force beyond the Debye screening length. 
The WKB approximation!^ provides an approximate 
treatment of the inhomogeneous nature of screening by 
using the local ionic strength; however, it does not cap¬ 
ture the long-range, accumulative effects. 

Since approximate treatments cannot fully account for 
all the features of inhomogeneous screening, previous cal¬ 
culations using these approximations for the double layer 
structure and interfacial properties are likely to be inac¬ 
curate. Such inaccuracy in treating the essential electro¬ 
static contributions makes it impossible to evaluate the 
relative importance of the various non-electrostatic ef¬ 
fects invoked, for example, to explain the surface tension 
behavior, such as the cavity energv [25 . hydration [2^, 
and dispersion forces 28|,l29|. In this Letter, we examine 
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the issue of inhomogeneous screening in salt solutions 
near a dielectric interface, by comparing the result from 
numerical solution of the full Green function with results 
obtained using approximate methods. We find that the 
effects of inhomogeneous screening on the double layer 
structure and interfacial properties are quite pronounced 
as the Debye screening length becomes comparable to the 
Bjerrum length. 

We have recently developed a general theory for weak- 
coupling systems with a fixed charge distribution pex in 
the presence of mobile cations with charge g+e and anions 
with charge q-e, in a dielectric medium of a spatially 
varying dielectric function e [lElii. The key result is the 
set of self-consistent equations for the mean electrostatic 
potential 'i/’C'') (nondimensionized by kT/e), the Green 
function G(r,r'), and the self-energy M±(r) of the ions: 


- V ■ (eVV^) = Pe. + rA+q+e-«+^-“+ - rA_g_e«-’^-“- 

( 1 ) 

- V • [eVG(r, r')] -b 2/(r)G(r, r') = d(r - r') (2) 


u±(r) = i J dr'dr"h±{r — r')G{r',r”)h±{r" — r) (3) 

where e = kTsoe/e^ is the scaled permittivity 
and A± is the fugacity of ions. T is introduced 
to constrain the mobile ions to the solvent region, 
/(r) = [( 7 ^c+(r)-b ( 7 ^c_(r)]/2 is the local ionic 

strength, with the ion concentration given by c±(r) = 
A±r exp [T9 ±V'(i’) “ '*^±(r)]. The short-range charge dis¬ 
tribution h±{r — r') on the ion is introduced to yield 
a finite Born solvation energy. Eq. [T] is the self-energy 
modified Poisson-Boltzmann equation, reflecting that the 
ion distribution is determined by both the mean electro¬ 
static potential and the self energy. From Eqs. [2] and 
[31 the inhomogeneity in the ionic strength affects the so¬ 
lution of the Green function and the self energy, which 
consequently affect the double layer structure through 
Eq. [TJ 

We now specify to a salt solution in contact with a low 
dielectric medium through a sharp interface (at z = 0) 
with fixed surface charge density Pex{r) = o’(-z)- Mobile 
ions are excluded from the low dielectric side. Both T 
and e are then step functions: T = 0 and e = ep for 
z < 0; r = 1 and e = es for z > 0. In the solvent region 
(z > 0), Eq. |T] becomes 

- = A+Q+e-«+^-“+ - A-g-e^-’^-^- (4) 

with boundary condition {dip/dz)z^o = —ajes- As¬ 
suming the solvent has a uniform dielectric constant in 
the entire z > 0 region, the Born energy is constant 
and can be absorbed into the reference chemical poten¬ 
tial. The remaining contribution is finite in the point- 
charge limit [h±{r — r') = q±S{r — r')], leading to the 


nontrivial and nonsingular part of the self energy as 
= (g±/2) lmv.^r [G(r,r') - 1/ (47res|r - r'|)]. 

To solve the Green function in the planar geometry, 
it is convenient to work in a cylindrical coordinate (r, z) 
and use the Fourier representation in the transverse di¬ 
rections: 


G(r,z,z') = 


27r 


kdkjQ{kr)G{k, z, z') (5) 


where Jq is the Oth-order Bessel function. It is easy to 
show that G(fc, z, z') satisfies: 


d G{k, 2 :, z ) p 2 / , 21 /^n i\ ^ r/' i\ 

- 5z2 -b [k (z) -b fc J G(fc, Z, z ) = —d{z, z ) 

( 6 ) 

for z > 0, with the boundary condition esdG/dz — 
kcpG = 0 at z = 0 37|. k(z) = [2/(z)/es]^^^ can be con¬ 
sidered the inverse of the local Debye screening length. 

The bulk screenin g a pproximation widely used in 
the literature [^, [UH^, replaces the spatially vary¬ 
ing screening length k{z) in Eq. [6] by the constant 
bulk screening length Kb, which enables an analyti¬ 
cal solution for the Green function as G{k,z,z') = 

/(^2esoj), where w = ^/Kf+¥ 

and A = (esw — epk)/{esu! + tpk). Substituting 
G{k, z, z') into Eq. pleads to the following intuitive form 
for the self energy when es ^ ep- 


<i± f , /e 

Ua- — T- I —Kb H- 


Sttes 


— 2KhZ 

Yz 


(7) 


with / = (e,9 — epl/fe? -j-ep) the dielectric contrast. The 
WKB approximation [39l441j is simply to replace the bulk 
Kb in Eq. |7|by the local «:(z) = [2/(z)/es]^^^. 

In this work, we perform the full numerical calcu¬ 
lation of the Green function using the finite difference 
method ^ |43[ . In order to ensure consistent numerical 
accuracy in removing the singularity of the same-point 
Green function, the free-space Green function satisfying 
—d^Gojdz^ + fc^Go = 5{z,z')/es, is also solved numeri¬ 
cally along with Eq. |6l is 


u±(z) = 


A 

47r 


G{k, z, z) — Go{k, z, z) kdk (8) 


Far away from the interface (z —> 00 ), the ion concentra¬ 
tion approaches the bulk value c^. It is straightforward 
to show A± = exp [—g^K{,/(87res)] [s^. 

We now apply the theory to salt solutions near the wa¬ 
ter/air interface (eg = 80 and ep = 1) with zero fixed 
surface charge (a = 0). This is the same system stud¬ 
ied by Onsager and Samaras. Figure 2(a) and 2(b) show 
the concentration profile of the ions in a 1:1 salt solution 
for two bulk concentrations. At low salt concentrations 
when ^ q^lp {Ip = l/dvreg is the Bjerrum length), 
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FIG. 2: Effect of inhomogeneous screening on the ion distribntion of (a)O.OlM and (b)l.OM 1:1 salt solution near the water/air 
interface, (c) the characteristic length of ion depletion vs. salt concentration. “Bulk k”, “WKB” and “Full Green” refer to the 
bulk screening approximation, the WKB approximation and numerically solving the full Green function, respectively. 


the effect of inhomogeneous screening is insignihcant, be¬ 
cause screening is weak even in the bulk. In this regime, 
both the bulk screening approximation and the WKB 
approximation are valid. 

In contrast, at higher salt concentrations when 
becomes comparable to or even smaller than inho¬ 
mogeneous screening affects the entire range of the de¬ 
pletion layer as shown in Fig. 2(b). Close to the interface 
{z < Ib), the ion concentration calculated by fully solv¬ 
ing the Green function is significantly lower than that 
predicted by the bulk screening approximation, because 
the local ionic strength that screens the image force is ob¬ 
viously smaller than the bulk. In Figs. 3(a) and 3(b) we 
provide a more visual representation of inhomogeneous 
screening by plotting the nondivergent part of the Green 
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FIG. 3: 2D visualization of the nondivergent Green function, 
G(r, r') — 1/ (47res|r — r'|), for l.OM 1:1 salt solution near the 
water/air interface. The test ion (white dot) is at z' = 1.5A 
for (a) and (b), and z' = 5A for (c) and (d). 


function, G(r,r') — 1/ (47res|r — r'|), i.e., the nondiver¬ 
gent part of the linear response electrostatic potential 
generated by point charge at a given distance from the 
interface. The potential generated by the ion close to 
the interface (z = 1.5A) is much stronger than that pre¬ 
dicted by the bulk screening approximation, the latter 
severely overestimating the local screening effect on the 
image charge interaction. Although this local effect is 
captured by the WKB approximation, neither of these 
two approximations capture the long-range and accumu¬ 
lative nature of the screening. The depletion layer cal¬ 
culated by fully solving the Green function extends to a 
range significantly longer than the bulk Debye screening 
length. As shown in Fig 3(d), even for an ion approach¬ 
ing the bulk solution {z = 5A, which is larger than the 
bulk screening length of 3.3A), the electric field from its 
image charge is not screened out. This remaining image 
charge interaction in turn has a long-range and accumu¬ 
lative effect that reinforces the field at the position of 
the point charge. The two approximate methods become 
progressively poorer as the salt concentration increases. 

We define d = [cb — c{z)\ dz/cb to characterize the 
width of the ion depletion layer, which is shown in Fig 
2(c) as a function of the salt concentration. The two 
approximate methods predict d to be an ever decreasing 
function of c^, determined by the bulk Debye screening 
length {d ~ In contrast, d calculated by fully solv¬ 

ing the Green function deviates significantly from the re¬ 
sults of the approximate methods as becomes compa¬ 
rable to q^ls, and reaches a constant value as Cb further 
increases up to IM. Thus, at high salt concentrations 
the image-charge repulsion renormalized by the inhomo¬ 
geneous screening creates a depletion layer of nearly con¬ 
stant width scaled by the Bjerrum length (d ~ q^ls) in 
stead of and becomes nearly independent of the salt 
concentration. 

As a consequence of the different behavior in the width 
of the depletion layer due to inhomogeneous screening, 
there is pronounced difference in the negative adsorption 
of ions (—F = [c{z) — Cb] dz) at the interface between 





















4 


results obtained by fully solving the Green function and 
those using the approximate methods as shown in Fig 
4. Because the approximate methods predict an ever de¬ 
creasing d as Cb increases , — F is a concave downwards 
function of Cb, as first shown by Onsager and Samaras [^. 
However, experimentally both —F and the surface ten¬ 
sion of the 1:1 salt solution increases essentially linearly 
with Cb in the range of O.IM < Cb < IM 21, To 
fit the experimental data, an ion-exclusion zone with 
constant width has been invoked in previous theoretical 
treatments 23|, [2^ . By fully accounting for the inhomo¬ 


geneous screening, our theory naturally predicts that —F 
increases linearly with Cb for O.IM < Cb < IM, as a con¬ 
sequence of a nearly constant d. In light of these results, 
it is quite possible that the inhomogeneous screening of 
the image force provides an explanation on the linear in¬ 
crease of the surface tension with the salt concentration; 
we are currently exploring this possibility. 



FIG. 4: Inhomogeneous screening effect on the negative ad¬ 
sorption of ions for a 1:1 salt solution at the water/air inter¬ 
face. 


Inhomogeneous screening has an even more pro¬ 
nounced effect on asymmetric salt solutions contain¬ 
ing multivalent ions, because multivalent ions are more 
strongly depleted than monovalent ions and are more ef¬ 
fective in screening. For a 2:1 salt solution, the divalent 
cations calculated by fully solving the Green function are 
pushed further away from the interface than predicted 
by the approximate methods as shown in Fig 5(a), lead¬ 
ing to a larger degree of charge separation. As a re¬ 
sult, the induced electrostatic potential is much larger 
than that obtained using the approximate methods; see 
Fig. 5(b). Such a large self-induced surface potential 
can significantly affect the interpretation of the zeta po¬ 
tential of colloidal surfaces^ and is a major contribu¬ 
tion to the Jones-Ray effect in the surface tension of salt 
solution 31, 4^4^. 

In conclusion, we have presented a self-consistent treat¬ 
ment of the inhomogeneous screening in salt solutions 
near a dielectric interface. The effect of inhomogeneous 



FIG. 5: Inhomogeneous screening effect on 0.05M 2:1 salt 
solution near the water/air interface, (a) Ion concentration 
scaled by c± and (b) dimensionless electrostatic potential. 


screening is twofold: the lower ionic strength near the 
interface results in less screening on the image force and 
hence stronger ion depletion, and the depletion zone has 
a long-range and accumulative effect on screening, which 
extends the range of the depletion layer. Gonsequently, 
the ion distribution is significantly affected when the bulk 
screening length is comparable to or smaller than the 
Bjerrum length. In this regime, the double layer struc¬ 
ture and the interfacial properties cannot be described 
by either the bulk screening approximation or the WKB 
approximation. The characteristic length of the deple¬ 
tion layer scales with the Bjerrum length, resulting in 
a linear increase of the negative adsorption of ions with 
concentration, in agreement with experiments. The in¬ 
homogeneous screening effect becomes is pronounced in 
the less polar solvent and for the ions of higher valency. 

Nonuniform ion distribution near a dielectric interface 
exists in many colloidal and biophysical systems. An ac¬ 
curate treatment of inhomogeneous screening is impor¬ 
tant to fully understand the role of electrostatic interac¬ 
tions, which in turn is necessary for evaluating the vari¬ 
ous nonelectrostatic contributions, such as the cavity en¬ 
ergy, hydration, and dispersion forces in these systems. 
The relative importance of these nonelectrostatic contri¬ 
butions, and even their existence, have been a subject of 
controversy [sij . By a more accurate treatment of the in¬ 
homogeneous screening effect, we are in a better position 
to understand some long-standing problems, such as the 
specific ion effects and salt concentration effects in the 
surface tension. 
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